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The semiclassical approximation has been used to study the effects of quantized fields on a 
classical background field in a wide variety of scenarios including black hole evaporation [T], the 
decay of an electric field due to the Schwinger effect [2] , heavy ion collisions in nuclear physics [3] , 
and preheating in chaotic inflation [H [5] . It is expected to be valid in cases where quantum effects 
are small, such as the initial stages of the evaporation of a solar mass black hole. At the opposite 
end of the spectrum is the case of preheating in models of chaotic inflation. Preheating occurs 
immediately after the inflationary phase and is a period in which the rate of particle production 
is extremely rapid, resulting in strong backreaction effects upon the inflaton field ilE]. It is not 
known whether the predictions of the semiclassical approximation can be trusted when quantum 
effects are so large. 

The semiclassical backreaction equations for quantum fields coupled to a classical background 
field arise out of the one loop effective action for that field [8]. As such they would typically be 
expected to break down when backreaction effects are large and terms coming from higher loops 
may be important. One way around this is to use a large N expansion where N is the number 
of identical quantum fields. The semiclassical backreaction equations become exact in the limit 
N ^ oo. This expansion has been used in cases such as preheating urn where backreaction effects 
are significant. 

For the semiclassical approximation to be valid, quantum fluctuations about the mean of what¬ 
ever quantity couples the quantum fields to the classical background field(s) must be small. One 
way to characterize these fluctuations in semiclassical gravity is through the two-point function for 
the energy-momentum tensor. However, for the symmetric part of this two-point function there 
can be state-dependent divergences [Tn|, and different renormalization schemes can yield different 
results when the points come together m To overcome these difficulties a criterion was given 
in |12j that relates the validity of the semiclassical approximation in gravity to the stability of so¬ 
lutions to the linear response equation, which results when the semiclassical backreaction equation 
is perturbed about a solution to that equation. The linear response equation has a term which 
involves the perturbed energy-momentum tensor, so renormalization proceeds in the usual way and 
there are no state dependent divergences. 

The criterion states that the large N semiclassical approximation in gravity will break down 
if any linearized gauge invariant quantity constructed from solutions to the linear response equa¬ 
tions with finite non-singular initial data, grows without bound. It has been shown to be satisfied 
for massive and massless free scalar fields in flat space in the Minkowski vacuum state [12], and 
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for conformally invariant free fields in the expanding part of de Sitter space, when spatially flat 
coordinates are used, the fields are in the Bunch-Davies state, and scalar perturbations are consid¬ 
ered m- Tensor perturbations for conformally invariant free fields were investigated in [H] and it 
was found that they are bounded, so the criterion is satisfied in that case as well. 

In this paper we continue an investigation begun in where we adapted the criterion in m to 
check the validity of the semiclassical approximation in models of preheating in chaotic inflation, in 
which rapid damping of the inflaton field occurs, and found evidence that quantum fluctuations are 
large in between the two periods of rapid damping. Here we study in great detail the relationship 
between solutions to the linear response equation and quantum fluctuations, and use the results 
to relate the size of the fluctuations to the particle production rate. We also include a case in 
which there is no rapid damping. We find evidence that quantum fluctuations are large and the 
semiclassical approximation breaks down whenever the particle production rate is high, including 
during the early stages of preheating when scattering effects ignored in our model and backreaction 
effects on the inflaton field are small. 

We consider a model of chaotic inflation for which the inflaton field cj) is coupled to N identical 
massless scalar fields 'i/’j with a coupling of the form ■ Full backreaction effects for this 

coupling have been investigated in detail in Refs. [U [Sj [E] (although not all of these were in the 
context of the large N expansion or for massless quantum fields). After a standard rescaling of the 
coupling constant g [5], the problem reduces to the coupling of the inflaton field to a single scalar 
field Tp, and the semiclassical backreaction equation for the inflaton field is — = 0. 

As in [5] we work in a flat space background and consider only homogeneous and isotropic solutions 
for (f). 

The mass of the inflaton field can be scaled out of the equations using t —>■ mt and </>—>■ p/m, 
with similar changes of variable for other relevant quantities. (See [5] for details.) The result is 

0 + (1 + = 0 ) (la) 





fk + + g'^(f^)fk = 0 . (ic) 

Note that is independent of the positive constant e and that M is a mass scale which typically 






4 


enters when computing renormalized quantities for massless fields [5]. 

The linear response equation can be derived as in [12] by taking a second variation of the effective 
action. The result is 


(□ - = 0 , 

= -ig'^ [ dV ct){x')5ct){x')e{t - t') {[i?{x),'il?{x')]) + ■ 


(2a) 

(2b) 


Here d(V’^)sD comes from a variation in the state of the quantum field. 

The linear response equation can also be derived by perturbing the semiclassical backreaction 
equation about one of its solutions. We illustrate this for homogeneous and isotropic perturbations. 


The equation for the inflaton field (la) and the mode equation (Ic) are perturbed in the usual way, 


keeping quantities that are first order in 5(j) and 6fk- The perturbed mode equation is then solved 
in terms of the solutions to ([^ with the result: 

Sfk = Akfk + Bkfk + 2gH f dt' [m)fk{t') - fk{t)mt')] ■ (3) 

Jo 

The coefficients Ak and are related to a change of state and are fixed by the initial values of 
6fk and its first derivative. Such a change in state (as pointed out in m for semiclassical gravity) 
must occur if the original state is a second order or higher adiabatic state |8]. 

The linear response equation in this case is 


6^ + il+ g"^+ g‘^(t)5{il?) = 0 , 

= 1^ [ dkk^ (fkSf*k + fk^fk) 


(4a) 


+ 


1 r 

47r^ 


dkk^ 


hsn+iisfk + J^'^ 




(4b) 


For the fourth order adiabatic states used in |5|, we find that = 0 to linear order. An explicit 
expression for Bj^ can easily be obtained but we will not display it here. 

If one can find solutions to ([^ then it is easy to generate approximate solutions to (2a). One 
simply takes two solutions, (pi and (j) 2 ^ which have nearly the same values at the initial time t = 0, 
and evolves them numerically in time. If we define the difference between the solutions to be 
5(pe ^ (t >2 — (pi: then 5(j)e satisfies the equation: 


64>e + (1 + g‘^{i>‘^)l)5(t)e + - (V'^)l)(01 + 5(pe 


(5) 
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The linear response equation (4a) in this case is 


54> + {l + g {ip )i)5(p + g 6{ip 4>i = 0 ■ 


( 6 ) 


Note that the first two terms in these equations have the same form. Thus, 5(p>e, which is a solution 
to ([^, is also an approximate solution to ([^ so long as the amplitude of the oscillations of cpi is 
much larger than the amplitude of oscillations of Scpe and S{ip‘^)\ « (V '^)2 — {'4’^)i- 

Because of the structure of Eq. ([^ it is possible to go further and separate out the part of the 
perturbation driven by 6{ip‘^). For simplicity choose the starting values for (p 2 and (pi such that 
</>2(0) = cp 20 , and ^ 2 ( 0 ) = (^i(O) = 0. Then let 


^<t>e = (p2- (pi = C(pi + 5(pc , 


(7) 


with c = {(p20 — (pio)/(pio- Substituting into (j^ and using (la) one finds that if 6(pe is an approximate 
solution to (©, then 6(pc is an approximate solution to the equation 


6(pc + {l + g‘^{ip‘^)i)d(pc + g‘^(pid{ip^)\ =-g‘^(pi6{ip‘^)\ . 

6(l)=S(f)c (50=c^i 


( 8 ) 


Thus the equation for 5(pc is the same as that for 6(p except that there is a source term which 
depends on S{ip'^) in (4b) evaluated with (p = (pi and 6(p = c(pi. Since the initial conditions are 


S(pc{0) = S(pc{0) = 0, at early times the growth is driven by the source term. 

As an illustration it is interesting to first look at a toy model in which {ip'^) in (la) is replaced 


by the last term in with e/M chosen so that this term is equal to —(p^j^ and the resulting 
equation for (p is (p + {1 — (p‘^)(p = 0. Then g‘^5{ip‘^) = —20^ and the source term for d(pc is 2c(p'^. 
The solutions for (p are stable for the starting values 0 < (p{0) < 1 and (p{0) = 0. In this case it is 
easy to solve the linear response equation directly. The results for (p{0) = 10“^ and = 10“® 

are shown in Fig. One sees that over the range shown there is linear growth in the amplitude of 
6(pc while the amplitude of 6(p does not grow significantly initially. This pattern of early growth of 
6(pc is also seen in the solutions to the full set of backreaction equations. 

In [6] it was predicted that there are two qualitatively different types of solutions to the back- 
reaction equation for (p. For one there is a relatively slow damping of the inflation field while for 
the other there is a period in which the inflaton field is rapidly damped. In [5] it was found for 
a flat space background that rapid damping of the inflaton field occurs whenever g‘^(pQ ~ 2 for 
models in which the starting values are (pQ = (p{t = 0) and (p{t = 0) = 0. Rapid damping does not 
occur for significantly smaller values such as g‘^(pQ = 1. Whenever rapid damping does occur it is 
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FIG. 1: Plotted are S(j) (upper curve) and S(j)c for the toy model described in the text with (/)(0) = 10 ^ and 
Scj)(0) = 10“®. The upper curve has been offset by 4 x 10“®. 

observed to happen twice and there appears to be no significant damping after that. These effects 
are illustrated in the upper panels of Fig. where the inflaton field is plotted as a function of time 
for = 1 and = 10. 

Examination of the plots in Fig. shows that in both cases 54>e grows exponentially at about 
the time that a significant amount of damping of cj) first occurs, while 5(j)c grows exponentially 
starting at much earlier times. After 5^c grows to be comparable in size to (5(/)e the two quantities 
are nearly identical and cannot be distinguished on the scale of the plots. For g‘^(l)Q = 1 a small 
amount of damping of (j) occurs very quickly followed by a much slower damping rate which goes on 
for a long time. During this latter period 54> grows approximately linearly in time. For g'^^'o = 10 
the exponential growth of 6(l)e continues through the end of the second rapid damping period and 
then all growth appears to cease. 

Using the detailed analysis of the particle production in [3] we find that the rate of growth of 
54>c appears to be closely tied to the overall particle production rate. It is exponential when the 
particle production rate is high, approximately linear in time during periods of slow damping when 
the rate is much smaller, and is negligible after the second rapid damping phase when the particle 
production rate is negligible (in cases where rapid damping occurs). 

In the case of preheating the growth rate of the solutions to the linear response equation varies 
significantly over time. Therefore the criterion in [12] should be modified so that the general form 
of the criterion is: in cases where a large N semiclassical approximation is used, the approximation 
will break down if any linearized gauge invariant quantity constructed from solutions to the linear 
response equations with finite non-singular initial data grows significantly for some period of time. 
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FIG. 2: The inflaton field </> is plotted versus time in the upper panels, and \5(j)e\ (dashed curves) and \5(j)c\ 
(solid curves) are plotted versus time in the lower panels. For each plot g = 10“^. For the upper left panel 
9^'Po = 1 s-iid tti® upper right one g‘^4>Q = 10. In the lower left panel 6(j)e is the difference between solutions 
with 3^00 = 1 + 10“^ and lower right panel 5(j)e is the difference between solutions with 

g^4>Q = 10(1 + 10“®) and = 10. 

For the model of preheating considered here, there are no gauge fields, so one can just consider 
solutions to the linear response equation. 

The quantity in ([^ may not be useful in all cases but it is useful in preheating where we have 
seen that it is a more sensitive measure of quantum fluctuations than 54> is. A second modification 
of the criterion specifically for preheating is that the semiclassical approximation breaks down if 
54>c grows rapidly for some period of time. It is clear from our results that this criterion is violated 
during the early stages of preheating, well before either scattering effects or backreaction effects 
are important. 

As pointed out in [6], the flat space approximation does not always give an accurate account 
of the details of the preheating process because the expansion of the universe can have a signifi- 








































































































































































































cant effect on the parametric amplification process. Nevertheless our resnlts strongly suggest that 
whenever there is a period in which a lot of parametric amplification occurs, the semiclassical 
approximation breaks down. 

This is the first application that has been made of the criterion in [12] for the validity of the 
semiclassical approximation when particle prodnction effects are significant. We think it likely, but 
cannot be certain, that our results generalize to similar situations and thus that the semiclassical 
approximation may never be valid when there is a high rate of particle production. 
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